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C$ • Abstract 



We construct an effective low-energy Hamiltonian from the clas- 
sical action via Monte Carlo with importance sampling. We use 
Monte Carlo (i) to compute matrix elements of the transition am- 
plitude and (ii) to construct stochastically a basis. The MC Hamil- 
tonian allows to obtain energies and wave functions of low-lying 
states. It allows also to compute thermodynamical observables in 
some temperature window (starting from temperature zero). We 
present examples from lattice field theory (Klein-Gordon model). 
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1 Introduction 

The most interesting phenomena in physics occur in many-body systems. 
The difficulty lies in solving models of such systems. A brilliant idea to do 
this is the renormalisation group approach a la Kadanoff, Wilson, Migdal 
and others. It suggests to compute critical phenomena by constructing 
an effective (renormalized) Hamiltonian, by thinning out degrees of free- 
dom in an appropriate manner. Here we suggest to adapt this idea for 
the pupose to construct an effective Hamiltonian in a low-energy window. 
The basic idea to thin out degrees of freedom is taken over from lattice 
field theory, where infinite-dimensional integrals (path integrals) are suc- 
cessfully simulated by a small number of representative configurations. In 
analogy to the representative configurations of path integrals, we construct 
a representative basis of Hilbert states via a Monte Carlo algorithm with 
importance sampling. Furthermore we use the Monte Carlo algorithm to 
compute matrix elements of the transition amplitude. 

In standard Lagrangian lattice field theory, one is used to construct 
a corresponding lattice Hamiltonian via the transfer matrix, which corre- 
sponds to the transition amplitude for infinitesimal short time (when the 
lattice spacing a t goes to zero). As a result, one finds a Hamilton operator 
(on the lattice) which still has the character of a many-body Hamiltonian. 
In contrast to that, we start out from the transition amplitude for a fi- 
nite time (in imaginary time corresponding to some finite temperature). 
We compute the transition matrix for all transitions, where the states are 
taken from the finite stochastically generated basis. As a result, we find 
an effective Hamiltonian, i.e. its spectrum and eigen states. It is evident 
that such Hamiltonian covers only a window of physics, which in this case 
is a temperature window, starting from temperature zero. 

From the viewpoint of the enormous success of standard Lagrangian 
lattice field theory, one may ask: Do we need such an effective Hamil- 
tonian? To give an example: Spectrum and wave functions of excited 
states. Wave functions in conjunction with the energy spectrum contain 
more physical information than the energy spectrum alone. Although lat- 
tice QCD simulations in the Lagrangian formulation give good estimates 
of the hadron masses, one is yet far from a comprehensive understand- 
ing of hadrons. Let us take as example a new type of hardrons made 
of gluons, the so-called glueballs. Lattice QCD calculations [[!]] predict 
the mass of the lightest glueball with quantum number J PC = ++ , to 
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be 1650 ± lOOMeV. Experimentally, there are at least two candidates: 
/o(1500) and /j(1710). The investigation of the glueball production and 
decays can certainly provide additional important information for experi- 
mental determination of a glueball. Therefore, it is important to be able to 
compute the glueball wave function. The standard Lagrangian lattice for- 
mulation has made slow progress also in some other areas, like on hadron 
structure functions (not only moments) in particular at small Q 2 and at 
small well as on scattering amplitudes. 



2 Effective Hamiltonian 



Let us discuss the construction of the effective Hamiltonian in several 
steps 0. First, consider in quantum mechanics in 1-dim the motion of a 
single particle of mass m under the influence of a local potential V(x). Its 
classical action is given by 

f fry 

S = J alt —x 2 - V(x) . (1) 

Given the classical action, one can determine Q.M. transition amplitudes. 
Like in Lagrangian lattice field theory, we use imaginary time in what 
follows. We consider the transition amplitude for some finite time T (tj n = 
0, tfi = T) and for all combinations of positions Xi,Xj. Here Xi, Xj run 
over a finite discrete set of points {xi, . . . , x^} located on the real axis. 
Suppose these points are equidistantly distributed (spacing Ax) over some 
interval [-L, +L]. The transition amplitudes are given by the (Euclidean) 
path integral, 



M^-(T) = [dx)exp[-S E [x]/h] 



Xi,T 

(2) 

Xj,0 

where Se denotes the Euclidean action. Suppose for the moment that 
those transition amplitudes were known. From that one can construct an 
approximate, i.e. effective Hamiltonian. Note that the matrix M(T) = 
[Mij(T)]NxN is a real, positive and Hermitian matrix. It can be factorized 
into a unitary matrix U and a real diagonal matrix D(T), 

M(T) = U ] D(T) U . (3) 

Because the matrix M(T) can be expressed in terms of the full Hamilto- 
nian H by 

M tJ (T) = ( Xi \e- HT / n \ Xj ) , (4) 



4 L.A. Caron, H. Kroger, G. Melkonyan, X.Q. Luo, K. Moriarty 



the matices U and D have the following meaning, 



U} k =< Xl \El ff > 
D k {T) = e- E l !!T ' h , 



(5) 



i.e., the k—th eigenvector \E1 > of the effective Hamiltonian H e ff can be 
identified with the k — th column of matrix W and the energy eigenvalues 
El of H e ff can be identified with the logarithm of the diagonal matrix 
elements of D(T). This yields an effective Hamiltonian, 



Note that in the above we have been mathematically a bit sloppy. The 
states \xi) are not Hilbert states. We have to replace \xi) by some "lo- 
calized" Hilbert state. This can be done by introducing box states. We 
associate to each Xi some box state 6j, defined by 



where Axi = Xi+i —X{. When we use an equidistant distribution of Xi, i.e. 
Axi = Ax, we refer to the basis of box states as the regular basis. Below 
we will consider a stochastic basis. 

3 Matrix elements 

We compute the matrix element My(T) directly from the action via Monte 
Carlo with importance sampling. This is done by writing each transition 
matrix element as a ratio of two path integrals. This is done by splitting 
the action into a non- interacting part and an interacting part, 



N 



Heff = E \K n > Et fS < E e k 



(6) 



k=i 




(7) 



S — Sq + Sy ■ 



(8) 



This allows to express the transition amplitude by 



M tJ (T)=M^(T) 



j+1 dz J[dx] exp[—Sv[x]/h] exp[—So[x]/h] 



(9) 
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Here e~ So ^ n is the weight factor and e~ Sv ^ h is the observable. M$ (T) 
stands for the transition amplitude of the noninteracting system, which 
is (almost) known analytically. For details see ref.[0. Carrying out these 
steps allows to construct an effective Hamiltonian, which reproduces well 
low energy physics, provided that the box functions cover an interval large 
enough (depending on the range of the potential) and the resolution Ax 
is small enough. This can be achieved with a small numerical effort for 
a one-body system in 1 dimension. Our goal is, however, to solve many- 
body systems. For such purpose, the above regular basis fails. What to 
do in such systems is the subject of the following section. 

4 Stochastic basis 

It is evident that the regular basis defined above becomes prohibitively 
large when applied to a many-body system. For example, in a spin model 
of a 1-dimensional chain of 30 atoms with spin 1/2, the Hilbert space 
has the dimension D = 2 30 = 1073741824. For such situations we wish to 
construct a smaller basis which gives an effective Hamiltonian reproducing 
well low-energy observables. Why should such a basis exist in the first 
place? The heuristic argument is the Euclidean path integral, which, when 
evaluated via Monte Carlo with importance sampling, gives a good answer 
for the transition amplitude. In particular, this is possible by taking into 
account a "small" number of configurations (e.g. in the order of 100 - 
1000). In a crude way the configurations correspond to basis functions. 
Thus we expect that suitably chosen basis functions exist, the number of 
which is in the order of 100 - 1000, which yields a satisfactory effective 
low energy Hamiltonian. 

How can we construct such a " small" basis? As example consider a free 
particle in D = 1 dimension. Recall: For the free system the transition 
amplitude reads 

G E ud(x } T; y, 0) = exp ^^T^ X ~ ^ ( 10 ) 

This is a positive function for all x, y, T. It can be interpreted as a 
probability density. We put y = and define a probability density P(x) 
by 

P{x) = ^G Eud {x,T;0,0), 
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Z = J dx G Eucl {x,T;0,0). (11) 



Then we define a selection process as follows: Using a random process 
with probability density P(x) one draws a "small" set of samples {x u \u e 
1, . . . , N e ff}. In the case of the free particle, P(x) is a Gaussian, 

PW ^ exp| -^' " = /f- (12) 

In other words, we select {x u } by drawing from a Gaussian distribution. 
We draw N e ff samples, giving N e ff states, where N e jf is considerably 
smaller than N, the original size of the basis. 

Next we ask: What do we do in the case when a local potential is 
present? The definition of the probability density P(x) given by Eq.(|TT|) 
generalizes to include also local potentials. In order to construct a stochas- 
tic basis one can proceed via the following steps: (i) Compute the Eu- 
clidean Green's function G E (x, t; 0, 0), e.g., by solving the diffusion equa- 
tion and compute P(x). (ii) Find an algorithm giving a random variable 
x distributed according to the probability density P(x) and draw sam- 
ples from this distribution, giving nodes, say x v . Finally, one obtains the 
stochastic basis by constructing the corresponding characteristic states 
from the nodes x v . The same goal can be achieved in an elegant and 
efficient manner via the Euclidean path integral. Writing Eq.flTTD as path 
integral yields 



x,T 

f[dy) exp[-S E [y]/h] 

P(x) = 



0.0 



xjT ■ (13) 

f+™dxf[dy)exp[-S E [y}/h] 



0.0 



Using a Monte Carlo algorithm with importance sampling (e.g., Metropo- 
lis) one generates representative paths, which all start at x = 0, t = 
and arrive at some position x at time t = T. Let us denote those paths 
(configurations) by C v = x u (t). We denote the endpoint of path C v at 
time t = T by x v = x u (T). Those form the stochastically selected nodes, 
which define the stochastic basis. 



5 Numerical results 
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5.1 Quantum mechanics 

The Monte Carlo Hamiltonian has been tested on a number of quantum 
mechanical potential models, by computing the spectrum, wave functions 
and thermodynamical observables in some temperature window |2|, |3|, |], |J. 
Although the results from the Monte Carlo Hamiltonian agree well with 
the exact results, low-dimensional models in Q.M. are not the target of 
this method. The target is rather high-dimensional models or models with 
a large number of degrees of freedom. The reason for this is based on the 
same principle which applies to numerical integration with Monte Carlo: 
The Monte Carlo method is not the most efficient method when doing low- 
dimensional integrals. However, it is most efficient when doing integrals 
in high dimensions (e.g. path integrals). 

5.2 Klein-Gordon model 

We consider in D=l a chain of coupled harmonic oscillators, which is 
equivalent to the Klein-Gordon field on a 1 + 1 lattice. Here we consider 
a chain of oscillators. The model is given by 

S = J dtT-V 

N 1 

T = Eo^n 

71=1 Z 
1 N 

V = ;,E^-^+i) 2 + «- (14) 

Z n=l 

The parameters have been chosen as m = 1, = 1, Qq = 2, N osc = 9, 
a = 1, T = 2 and h = 1. For the adjustable parameter a in the stochastic 
basis, we choose a = 'ftsinh(fiT) / (mfi). After the stochastic basis with 
N s toch — 1000 is generated, we obtain the matrix elements M n > n Then we 
compute the eigenvalues and eigenvectors using the method described in 
Sect. Table [l] gives a comparison between the spectrum from the effective 
Hamiltonian with the stochastic basis and the analytic result for the first 
20 states. They agree very well. This means that the Gaussian distribu- 
tion for the stochastic basis works well. 

We have also computed thermodynamical quantities such as the partition 
function Z, free energy F, average energy U = E and specific heat C. 
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Table 1: Spectrum of Klein-Gordon model on the lattice, MC Hamiltonian 
(stochastic basis) versus exact result. 



n 




/ricxact 


1 


10.904663192168 


10.944060480668 


2 


12.956830557334 


12.944060480668 


3 


12.985023578737 


13.057803869484 


4 


13.044311582647 


13.057803869484 


5 


13.299967341242 


13.321601993380 


6 


13.345480638394 


13.321601993380 


7 


13.552195133687 


13.589811791733 


8 


13.585794986361 


13.589811791733 


9 


13.680136748933 


13.751084748745 


10 


13.744919087477 


13.751084748745 


11 


14.984737011385 


14.944060480668 


12 


15.012353803145 


15.057803869484 


13 


15.057295761044 


15.057803869484 


14 


15.108904652020 


15.171547258300 


15 


15.125356713561 


15.171547258300 


16 


15.187413290039 


15.171547258300 


17 


15.308536490102 


15.321601993380 


18 


15.396255686587 


15.321601993380 


19 


15.420708031412 


15.435345382196 


20 


15.432823810789 


15.435345382196 



The analytical results are 

iVosc 1 

Z W = Tt(exp M= n 2sinh(M . 

E{(3) = -Tr( J r7exp(-/3 J f7)) = -^-=5:^coth(^/2), 

n,o^ , 9E , o2 dE , { phui/2 \ 2 

°^ = k °W = ~ k ^ W = g ( 2^(^/2) J ' (15) 



where 



Ui = ^n 2 + 4n 2 sin 2 ( Pl Ax/2) , 
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Ap = 2n/(N osc Ax) , 
Xj = [-(N osc -l)/2 + (j-l)]Ax , 
Pi = [-(N osc -l)/2 + (l-l)]Ap. (16) 

Here j and I run from 1 to iV osc (number of oscillators). Ax = a = 1 is the 
lattice spacing, (5 = T/h, the temperature is related to j3 via T = l/([3k B ), 
and ks is the Boltzmann constant. Since we have approximated H by H c s, 
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Figure 1: Specific heat (C/ks) of the Klein-Gordon model on a 1+1 dimensional 
lattice. 

we can express those thermodynamical observables via the eigenvalues of 
the effective Hamiltonian 

N 



n=l 



E (;3) = tr^w 

( N ( E eS\2 -pEf 

C7-C9) = ksffc^*^ (F" 



\n=l 



(17) 



Since this is a static system, the eigenvalues should in principle not vary 
with p. This is observed numerically within statistical errors when j3 and 
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N stoch are not too small. Using this assumption and the spectrum at 
f3 — 2, we can compute thermodynamical quantities for other values of f3. 
The specific heat as a function of f3 is shown in Fig. [1]. Again, the results 
from the MC Hamiltonian are in good agreement with the analytical ones 
when /3 > 1. Preliminary results for larger 7V osc indicate that it is not 
necessary to increase N stoc h accordingly. This property is very important 
for an application of the algorithm to many-body systems and QFT. 

6 Summary 

In this paper, we have extended the effective Hamiltonian method with 
a stochastic basis to QFT, and taken the Klein-Gordon model as an ex- 
ample. The results are very encouraging. We believe that the application 
of the algorithm to more complicated systems will be very interesting. It 
will allow a non-perturbative investigation of physics beyond the ground 
state. 
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